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Many inflating spacetimes are likely to violate the weak energy condition, a key assumption of 
singularity theorems. Here we offer a simple kinematical argument, requiring no energy condition, 
that a cosmological model which is inflating — or just expanding sufficiently fast — must be incomplete 
in null and timelike past directions. Specifically, we obtain a bound on the integral of the Hubble 
parameter over a past-directed timelike or null geodesic. Thus inflationary models require physics 
other than inflation to describe the past boundary of the inflating region of spacetime. 


PACS numbers: 98.80.Cq, 04.20.Dw 


I. Introduction. Inflationary cosmological models fil gg 
are generically eternal to the future rl . In these mod- 
els, the Universe consists of post-inflationary, thermal- 
ized regions coexisting with still-inflating ones. In co- 
moving coordinates the thermalized regions grow in time 
and are joined by new thermalized regions, so the comov- 
ing volume of the inflating regions vanishes as t — oo. 
Nonetheless, the inflating regions expand so fast that 
their physical volume grows exponentially with time. As 
a result, there is never a time when the Universe is com- 
pletely thermalized. In such spacetimes, it is natural 
to ask if the Universe could also be past-eternal. If it 
could, eternal inflation would provide a viable model of 
the Universe with no initial singularity. The Universe 
would never come into existence. It would simply exist. 


This possibility was discussed_in the early days of in- 
flation, but it was soon realized that the idea could 
not be implemented in the simplest model in which the 
inflating universe is described by an exact de Sitter space. 
More general theorems showing that inflationary space- 
times are geodesically incomplete to the past were then 
proved One of the key assumptions made in these 
theorems is that the energy-momentum tensor obeys the 
weak energy condition. Although this condition is satis- 
fied by all known forms of classical matter, subsequent 
work has shown that it is likely to be violated by quan- 
tum effects in inflationary models (9, fia}. Such viola- 
tions must occur whenever quantum fluctuations result 
in an increase of the Hubble parameter H — i.e., when 
dH/dt > 0 — provided that the spacetime and the fluc- 
tuation can be approximated as locally flat. Such upward 
fluctuations in H are essential for the future-eternal na- 
ture of chaotic inflation. Thus, the weak energy condi- 
tion is generally violated in an eternally inflating uni- 
verse. These violations appear to open the door again to 
the possibility that inflation, by itself, can eliminate the 
need for an initial singularity. Here we argue that this is 
not the case. In fact, we show that the general situation 


is very similar to that in de Sitter space. 

The intuitive reason why de Sitter inflation cannot be 
past-eternal is that, in the full de Sitter space, expo- 
nential expansion is preceded by exponential contraction. 
Such a contracting phase is not part of standard inflation- 
ary models, and does not appear to be consistent with 
the physics of inflation. If thermalized regions were able 
to form all the way to past infinity in the contracting 
spacetime, the whole universe would have been thermal- 
ized before inflationary expansion could begin. In our 
analysis we will exclude the possibility of such a con- 
tracting phase by considering spacetimes for which the 
past region obeys an averaged expansion condition, by 
which we mean that the average expansion rate in the 
past is greater than zero: 


Hay > 0. (1) 


With a suitable definition of H and the region over which 
the average is to be taken, we will show that the averaged 
expansion condition implies past-incompleteness. 

It is important to realize that the terms expansion and 
contraction refer to the behavior of congruences of time- 
like geodesics (the potential trajectories of test particles). 
It is meaningless to say that a spacetime is expanding at 
a single point, since in the vicinity of any point one can 
always construct congruences that expand or contract at 
any desired rate. We will see, however, that nontrivial 
consequences can result if we assume the existence of a 
single congruence with a positive average expansion rate 
throughout some specified region. 

While the past of an inflationary model is a matter of 
speculation, the attractor nature of the inflationary equa- 
tions implies that many properties of the future can be 
deduced unambiguously. According to the standard pic- 
ture of inflation, all physical quantities are slowly vary- 
ing on the scale of H~!. In the vicinity of any point P 
in the inflating region, we can choose an approximately 
homogenous, isotropic and flat spacelike surface which 


can serve as the starting point for the standard analy- 
sis of stochastic evolution fi}. A simple pattern of ex- 
pansion is established, in which the comoving geodesics 
x = const in the synchronous gauge form a congruence 
with H > \/(87/3)Gpo, where po is the minimum energy 
density of the inflationary part of the potential energy 
function. This congruence covers the future light cone 
of P. While large fluctuations can drive H to negative 
values, such fluctuations are extremely rare. Once in- 
flation ends in any given region, however, many of the 
geodesics are likely to develop caustics as the matter 
clumps to form galaxies and black holes. If we try to 
describe inflation that is eternal into the past, it would 
seem reasonable to assume that the past of P is like the 
inflating region to the future, which would mean that a 
congruence that is expanding everywhere, except for rare 
fluctuations, can be defined throughout that past. 

For the proof of our theorem, however, we find that it is 
sufficient to adopt a much weaker assumption, requiring 
only that a congruence with H,, > 0 can be continu- 
ously defined along some past-directed timelike or null 
geodesic. 

In Section II, we illustrate our result by showing how 
it arises in the case of a homogeneous, isotropic, and spa- 
tially flat universe. In the course of the argument we shed 
some light on the meaning of an incomplete null geodesic 
by relating the affine parameterization to the cosmologi- 
cal redshift. In Section III we present our main, model- 
independent argument. In Section IV we offer some re- 
marks on the interpretation and possible extensions of 
our result. 


II. A simple model. Consider a model in which the metric 
takes the form 


ds? = dt? — a?(t)dz#?. (2) 


We will first examine the behavior of null geodesics, and 
then timelike ones. 

From the geodesic equation one finds that a null 
geodesic in the metric (2), with affine parameter À, obeys 
the relation 


dX x a(t) dt. (3) 


Alternatively, we can understand this equation by consid- 
ering a physical wave propagating along the null geodesic. 
In the short wavelength limit the wave vector k” is tan- 
gential to the geodesic, and is related to the affine pa- 
rameterization of the geodesic by k” œ dz” /dà. This al- 
lows us to write dA x dt/w, where w = k? is the physical 
frequency as measured by a comoving observer. In an ex- 
panding model the frequency is red-shifted as w « 1/a(t), 
so we recover the result of Eq. (8). 

From Eq. B), one sees that if a(t) decreases suffi- 
ciently quickly in the past direction, then f a(t) dt can be 
bounded and the maximum affine length must be finite. 


To relate this possibility to the behavior of the Hubble 
parameter H, we first normalize the affine parameter by 
choosing dà = [a(t)/a(ty)] dt, so dA/dt = 1 when t = ty, 
where ty is some chosen reference time. Using H = å/a, 
where a dot denotes a derivative with respect to t, we can 
multiply Eq. B) by H(A) and then integrate from some 
initial time t; to the reference time tf: 


Alts) alts) da 
f H(A)dà = / zi, (4) 
A(ti) a(ti) alts) 


where equality holds if a(t;) = 0. Defining Hav to be an 
average over the affine parameter, 


1 ts Bhai ee 1 
wont, HOS say 
(5) 
we see that any backward-going null geodesic with Ha, > 
0 must have a finite affine length, i.e., is past-incomplete. 
A similar argument can be made for timelike geodesics, 
parameterized by the proper time 7. For a particle of 
mass m, the four-momentum P” = m dx” /dr, so we can 
write dr = (m/E)dt, where E = P° is the energy of 
the particle as measured by a comoving observer. If we 
define the magnitude of the three-momentum p by p? = 
—gi; P! PÍ, where i and j are summed 1 to 3, then E = 
Vp? +m?. For a comoving trajectory we have Pt = 0, 
and therefore dr = dt. For all others, p x 1/a(t) fd, so 
we can write p(t) = [a(t;)/a(t)] pf, where pr denotes the 
value of the three-momentum p at the reference time tz. 
Combining all this, we find 


Haw 


Il lll 
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where the inequality becomes an equality when a(t;) = 0. 


Here Ey = \/ pe +m? and y = 1/,/1—v22,, where 


Urel = pr/Eyz is the speed of the geodesic relative to 
the comoving observers at time tr. Since the integral 
is bounded, the argument used for null geodesics can be 
repeated, with the average taken over proper time. 


II. The main argument. In this section we show that 
the inequalities of Eqs. (£) and (@) can be established in 
arbitrary cosmological models, making no assumptions 
about homogeneity, isotropy, or energy conditions. To 
achieve this generality, we need a definition of the Hubble 
parameter H that applies to arbitrary models, and which 
reduces to the standard one (H = a/a) in simple models. 
Consider a timelike or null geodesic O (“the observer” ). 
We assume that a congruence of timelike geodesics (“co- 
moving test particles” ) has been defined along O ig, and 
we will construct a definition for H that depends only on 
the relative motion of the observer and test particles. 


Test Test 
Particle 2 


Particle 1 


FIG. 1: The observer’s worldline O and two test particles. 


In order to motivate what we do, we first consider 
the case of nonrelativistic velocities in Minkowski space. 
Suppose that the observer measures the velocities of the 
test particles as a function of the time 7 on his own clock. 
At time 7; particle 1 passes with velocity w(71), and at 
time T2 = 7, + Ar particle 2 passes with velocity (r2). 
What expansion rate could he infer from these measure- 
ments? The separation vector between the positions r 
and fə of the two particles at T2 is AF =r, — 7 = U AT, 
and its magnitude is Ar = |Ar'|. Their relative velocity 
is AU = t — üz = —(dti/dr)Ar. The Hubble expansion 
rate is defined in terms of the rate of separation of these 
particles, which in turn depends on the radial compo- 
nent of their relative velocity, Au, = At - Ar/Ar. The 
inferred Hubble parameter H is then 


Au, Aŭ. Ar 
= Dr aFe Pe? 


which will equal the standardly defined Hubble parame- 
ter for the case of a homogeneous, isotropic universe. The 
expression for H may be simplified to H = —d In vye1/dr, 
where vrei = |u |. The fact that H is the total derivative 
of a function of vre} implies that the variation of Ure] is 
determined completely by the local value of H, even if 
the universe is inhomogeneous and anisotropic. 

We can now generalize this idea to the case of arbitrary 
velocities in curved spacetime. Let v“ = dx“ /dr be the 
four-velocity of the geodesic ©, where we take 7 to be 
the proper time in the case of a timelike observer, or an 
affine parameter in the case of a null observer. Let u/(r) 
denote the four-velocity of the comoving test particle that 
passes the observer at time 7. We define y = u,v”, so in 
the timelike case y may be viewed as the relative Lorentz 
factor (1/4/1 — v2.) between u” and v”. In the null case, 
y = dt/dr, where t is the time as measured by comoving 
observers, and 7 is the affine parameter of O. 

Consider observations made by © at times 7, and 
T = 7, + Ar, as shown in Fig. [] where Ar is infinites- 
imal. Let Ar“ be a vector that joins the worldlines of 


the two test particles at equal times in their own rest 
frame. Such a vector is perpendicular to the worldlines, 
and can be constructed by projecting the vector —v“Ar 
to be perpendicular to u”: Art = —v!Ar + yu”Ar. 
In the rest frame of the observer this answer reduces 
for small velocities to its nonrelativistic counterpart, 
Ar! = (0,u' Ar). This is a spacelike vector of length 
Ar = |Ar"| = y7? — K Ar, where k = vv, is equal 
to 1 for the timelike case and 0 for the null case. The 
separation velocity will be Au” = —(Du*/dr) Ar, where 
D/dr is the covariant derivative along O. The covariant 
derivative allows us to compare via parallel transport vec- 
tors defined at two different points along O, and can be 
justified by considering the problem in the free-falling 
frame, for which the affine connection vanishes on O at 
T = Tı. The radial component of this velocity will be 
Au, = — (Au"Ar,,) /Ar, where the sign arises from the 
Lorentz metric. We define the Hubble parameter as 


H= Au, — —v,(Du"/dr) 


i wee (8) 


Since O is a geodesic, we have (Dv“/dr) = 0, and 
therefore 


—dy/dr d 
y= UA Ap 
-k dr (0)), (9) 
where 
yo} null observer (x = 0) 
By) = 44, 1 (10) 
z! oe timelike observer (K = 1) 
y— 


As in Section II, we now integrate H along O from some 
initial 7; to some chosen T+: 


f ar= FoF) Fo). (a) 


In the null case F (yf) = Yes which is equal to the value 
of dr/dt at tf, normalized in Section II to unity. 

Eq. (I) therefore reproduces exactly the results of 
Eqs. (£) and (6), but in a much more general context. 
Again we see that if H,, > 0 along any null or non- 
comoving timelike geodesic, then the geodesic is neces- 
sarily past-incomplete. 


IV. Discussion. Our argument shows that null and time- 
like geodesics are, in general, past-incomplete in infla- 
tionary models, whether or not energy conditions hold, 
provided only that the averaged expansion condition 
Hav > 0 holds along these past-directed geodesics. This 
is a stronger conclusion than the one arrived at in pre- 
vious work in that we have shown under reasonable 
assumptions that almost all causal geodesics, when ex- 
tended to the past of an arbitrary point, reach the bound- 
ary of the inflating region of spacetime in a finite proper 
time (finite affine length, in the null case). 


What can lie beyond this boundary? Several possibil- 
ities have been discussed, one being that the boundary 
of the inflating region corresponds to the beginning of 
the Universe in a quantum nucleation event . The 
boundary is then a closed spacelike hypersurface which 
can be determined from the appropriate instanton. 

Whatever the possibilities for the boundary, it is clear 
that unless the averaged expansion condition can some- 
how be avoided for all past-directed geodesics, inflation 
alone is not sufficient to provide a complete description of 
the Universe, and some new physics is necessary in order 
to determine the correct conditions at the boundary Lg. 
This is the chief result of our paper. The result depends 
on just one assumption: the Hubble parameter H has a 
positive value when averaged over the affine parameter 
of a past-directed null or noncomoving timelike geodesic. 

The class of cosmologies satisfying this assumption is 
not limited to inflating universes. Of particular interest is 
the recycling scenario LA, in which each comoving region 
goes through a succession of inflationary and thermalized 
epochs. Since this scenario requires a positive true vac- 
uum energy Py, the expansion rate will be bounded by 
Hmin = y 8TGp/3 for locally flat or open equal-time 
slicings, and the conditions of our theorem may be satis- 
fied. One must look carefully, however, at the possibility 
of discontinuities where the inflationary and thermalized 
regions meet. This issue requires further analysis. 

Our argument can be straightforwardly extended to 
cosmology in higher dimensions. For example, in the 
model of Ref. brane worlds are created in collisions 
of bubbles nucleating in an inflating higher-dimensional 
bulk spacetime. Our analysis implies that the inflating 
bulk cannot be past-complete. 

We finally comment on the cyclic universe model 
in which a bulk of 4 spatial dimensions is sandwiched 
between two 3-dimensional branes. The effective (3 + 1)- 
dimensional geometry describes a periodically expanding 
and recollapsing universe, with curvature singularities 
separating each cycle. The internal brane spacetimes, 
however, are nonsingular, and this is the basis for the 
claim that the cyclic scenario does not require any 
initial conditions. We disagree with this claim. 

In some versions of the cyclic model the brane space- 
times are everywhere expanding, so our theorem imme- 
diately implies the existence of a past boundary at which 
boundary conditions must be imposed. In other versions, 
there are brief periods of contraction, but the net result of 
each cycle is an expansion. For null geodesics each cycle 
is identical to the others, except for the overall normal- 
ization of the affine parameter. Thus, as long as Hav > 0 
for a null geodesic when averaged over one cycle, then 
Hav > 0 for any number of cycles, and our theorem would 
imply that the geodesic is incomplete. 
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